Probing Low-x QCD With Very High Energy Prompt Muons 
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We explore the possibility of utilizing the prompt muon fluxes at very high energies in order to 
discriminate various models/parametrizations of low-x QCD behaviour of hadronic cross-sections 
relevant at such energies. We find that the pair meter technique for measuring high energy prompt 
muons can be very efficient in such an endeavor. As a by product, it allows to cleanly probe the 
change in composition of the primary cosmic rays expected at high energies. 
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Introduction The cosmic ray (CR) spectrum is 
characterized by a sharply falling power law behaviour, 
^ ~ £; _ (t+ 1 ) [1]. The spectrum gets more steeper 
around 10 6 GeV with the spectral index 7 changing 
from 1.7 to 2.1 - this region is called the knee. Around 
£~5x 10 9 GeV, one observes a flattening of the spec- 
trum, with the spectral index 7 falling between 1.4 and 
1.7. This is the so called ankle. The change in the slope 
of the spectrum at these two places is a puzzling issue. 
The region beyond the ankle is the regime of ultra high 
energy cosmic rays. There is not much data available in 
that region and no clear consensus exists on the com- 
position or the particle content in this region [3]. It is 
generally believed that the change in the slope around 
the knee is astrophysical in nature rather than any spe- 
cific change in hadronic properties and/or interactions. 
Though not conclusive, there appears to be some evi- 
dence [3( that the composition is heavier above the knee 
region. If true, then a significant suppression of the very 
high energy (> 10 5 GeV) lepton fluxes is expected due 
to CR interactions during the journey downwards [1]. 

There is a sharp reduction of the lepton flux from pi- 
ons and kaons (this component is called the conventional 
flux) above a few TeV [5[ . This is due to the increasing 
competition between their interaction and decay lengths. 
Therefore, at very high energies, almost all the lepton 
flux is expected to arise from the semi-leptonic decays 
of heavy hadrons (this component is called the prompt 
flux), most noticeably the charmed hadrons. It must 
be mentioned that the B-hadrons can also contribute 
significantly to the lepton fluxes in this energy range. 
This perhaps becomes an important issue for the v T be- 
cause the only source of v T in the charm sector is D s 
while there can be many B-hadron decay channels giv- 
ing rise to r and v T . Another point of importance is the 
fact that the neutrinos propagating at such high energies 
have almost vanishing probabilities to oscillate to an- 
other flavour. Therefore, measurement of neutrino flux 
in this regime requires no corrections due to the oscilla- 
tion phenomenon. Further, it is known that the prompt 
muon flux is only about 10% smaller than the prompt 
Vfj, flux at the surface of the earth. Therefore, measure- 
ment of the prompt muon flux at high energies will act 
as a normalization for the prompt neutrino flux and a 
direct comparison of the two is both desirable and also 
necessary. The prompt muon flux, thus, is an object 
of great importance, not just from the above mentioned 



points but also from the fact that such muons are a ma- 
jor background to various neutrino experiments. 

To this end, we require precise theoretical predictions 
for the prompt lepton fluxes. However, in reality the sit- 
uation is drastically the opposite. Various phenomcno- 
logical predictions for the prompt lepton fluxes span over 
two decades of interval [1 E H, B EE [H O • The reason 
is mainly due to the vastly different choices for the charm 
production cross-section - perturbative QCD (pQCD) 
with a K factor next-to-leading order (NLO) pQCD 
flit - quark-gluon string models and recombination 
quark-parton models [l(| • Within the QCD based mod- 
els, the problem can be traced back to the fact that the 
available data from accelerator/collider experiments is 
used as an input for the charm production (and also 
for other hadrons) at those experiments. This is then 
extrapolated to high energies and low-x values. Typical 
values sampled in high energy cosmic ray events, relevant 
for prompt muon flux via the charm production, corre- 
spond to x ~ 1CP 5 and maybe smaller. There is no data 
available in this region and one is forced to extrapolate 
without clear guidelines. The gluon distribution func- 
tion g(x) extracted from the available data, and naively 
extrapolated to higher energies and lower x values, is 
known to grow indefinitely, thereby causing concern re- 
garding unitarity. Physically, it is expected that some 
mechanism should tame this sharply rising behaviour, 
though none is conclusively known at present. Often, 
the theoretical models assume a power law behaviour 
for the gluon distribution function, xg(x) ~ x~ x , with 
A varying in the range — 0.5. The fluxes therefore 
depend strongly on the chosen value of A. The data 
from ep scattering at HERA, albeit at a very different 
energy and low-x regime, show a rise in the gluon dis- 
tribution function. On the other hand, RHIC data at 
higher energies and for different compositions of the col- 
liding particles (also later HERA data) gives some indi- 
cation of partial taming of the rising distribution func- 
tion. This is the so called saturation of gluon density 
at very low-x. Roughly speaking, at very low-x values, 
one can expect that the non-linear effects, recombina- 
tion effects to become significant, thereby leading to the 
saturation or taming of the gluon density. For an early 
discussion of RHIC results and theoretical expectations, 
see [l3j]. The saturation models have been shown to fit 
the HERA data as well. An important lesson is to be 
very cautious in extrapolating the existing accelerator 
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data to cosmic ray energies. The saturation picture can 
be expected to be applicable at cosmic ray energies fl4j. 
For a general overview of various theoretical issues in 
low-x and saturation physics see [HI]. An overview of 
hadronic interactions and cosmic rays can be found in 

The experimental situation is not very precise either 
at this stage. Various experiments [17| provide upper 
limits on the muon and neutrino fluxes in the energy 
range of interest. At present, these limits allow a large 
variation in the prompt fluxes. One can therefore ex- 
pect that better measurements of muon fluxes can play 
a definitive role in selecting the charm production mod- 
els, and thereby, also providing invaluable information 
about parton densities at such low-x and high energy 
values. Another related source of large theoretical un- 
certainties is strong dependence of the hadronic cross- 
sections on the renormalization and factorization scales. 
This is partly related to the naive extrapolation of par- 
ton distribution functions to very different energy and 
x-values. For the case of conventional fluxes originating 
from the pions and kaons, these issues are in much better 
control and therefore the predictions stand on a sound 
footing. 

In the present letter, we explore the possibility of uti- 
lizing the high energy prompt muon flux(es) in order to 
investigate whether the general expectations expressed 
above can in practice help in selecting the charm pro- 
duction model/parametrization. The range of models 
or parameterizations employed in the literature is too 
vast to be covered in this work. We choose some of 
the models often used and compare the predictions, in- 
corporating the saturation model of Golec-Biernat and 
Wuthsoff [18(]. Also, it should be mentioned that for the 
case of muons and , the flux from the beauty hadrons 
is not more than a few percent of that from the charm 
hadrons. Therefore, we ignore the beauty contributions 
in this work. However, for v T flux, there is almost 40% 
enhancement [l2| and should be considered. 

Prompt muons and their detection At higher 
and higher energies, the (prompt) muon flux is expected 
to reflect the onset of contributions from the production 
of heavy hadrons - the charmed and the beauty. The en- 
ergy dependence of prompt muons (the flux is isotropic) 
follows the parent cosmic ray spectrum and is therefore 
harder. In contrast, the conventional flux energy spec- 
trum is more steeply falling (almost by one extra power), 
caused due to the tension between the decay and inter- 
action of the parent mesons while passing through the 
atmosphere. A schematic picture of the entire process, 
very similar to the conventional lepton flux, is as follows: 
CR flux(X = 0) — > Flux(X) — > Charm production — > 
— > Decay — > Prompt Flux — > To detector 
Of all the ingredients entering the calculation, the model 
for charm production has the biggest uncertainties as dis- 
cussed above. Before proceeding further, let us have a 
look at the expression for the lepton flux. Assuming in- 
finite isothermal atmospheric depth, the final expression 
for the flux, for lepton energies Ei < e-^iarm K ^ *-*eV 
is (for zenith angle < 60°, ef u = 7 -^, h = 6.4 km 
where nij , Ti are the rest mass and life time of the ith 



particle) 



Zyj(E) Z NYc (E) ^4It 4n(E,0) 



Y a ^l N-\-A — >-Y c 



Xn{E^ 

depletion 



(i) 



where A N (E) = , X ^ {E ] P , and X N (E) = patm( *~\ 

« v j 1-Z NN (E) " \ I a NA n A (X=h) 

are the attenuation and interaction length of the nucleon 
in the atmosphere. <Pn(E,0) is flux of the nucleon at 
depth X=0. Y c denotes various charm hadrons in the 
sequence - D ^ , D s , A c . For the total Nucleon- Air cross 
section, we adopt the parametrization [1 91 ] 



o~na(E) = 



280 - 8.71n(^) + 1.14 In 2 (^) 



(2) 

Assuming scaling, spectrum weighted moments for a 
power law initial flux, <j>N{E,X = 0) = (/)onE n ^ +1 \ 
are defined as 



Z Nj (E) = f 
Jo 



1 



dxx 1 1 
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(3) 



where x ~ -4r and j denotes any of the final species. 
The charm hadron Z moments are related to the par- 
tonic level charm Z moments through their relative yields 
EE?], Z Nj (E) = fjZ Nc (E). The lepton production Z 
moments can be defined in an analogous fashion by re- 
placing the production cross-section by the (differential) 
decay rates and plugging in the corresponding branching 
fractions [2l| . 

It is therefore quite clear from the above expressions 
that the lepton fluxes at the end are quite sensitive to 
the charm production cross section. Till the knee, the 
cosmic ray flux and composition is rather established 
and therefore, the only source of large error is the charm 
cross section. This therefore gives us a unique possibil- 
ity to gain information about heavy quark production 
mechanism at high energies and low-x. Due to very pen- 
etrating nature of the muons, it is extremely hard to have 
a direct detection. The direct measurement of such high 
energy muons will require impractically large detectors 
with, presumably, very strong magnetic fields. Also the 
error in magnetic deflection of the muons at such large 
energies become comparable to the total angular deflec- 
tion itself. However, there exists a clean method capable 
of capturing the energy information of the incoming high 
energy muons. The method, pair meter technique [221 ] . is 
based on measuring the individual muon energy via the 
energy deposition by fast muons in the dense matter. It 
is quite well known that above a few TeVs, the dominant 
energy loss mechanism for muons while passing through 
matter is the e + e~ pair production, valid for the whole 
range of momentum transfer values. Further, at TeV 
energy or above, the differential energy loss, dE^/dx is 
linear in E^, providing a clean estimation of muon en- 
ergy. This technique has been successfully employed for 
a limited cosmic ray data [23| and the study confirms 
that the same method is capable of being used for even 
high energy muons. The requirement for such a mea- 
surement is to have the muons pass through a dense ma- 
terial like iron (or some other heavier material). In this 
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regard, it is worthwhile to emphasize that the typical ex- 
pressions employed in most cosmic ray studies are Born 
cross sections. However, for heavier materials, there is a 
need to consider the Coulomb corrections arising due to 
all order resummation effects coming from Za em , where 
Z is the atomic number of the material [24(. One can 
easily verify that for the case of iron, even at such large 
muon energies and for small momentum transfer, these 
corrections are at most a couple of percent. For lead, on 
the contrary, the Coulomb corrections are way too large 
and can not be neglected. Therefore, iron seems to be 
ideal for such a case. We, thus, choose to investigate the 
capability of this method in context of an iron calorime- 
ter. For the sake of illustration we choose to work with 
a volume of 50 kT. This choice is not completely arbi- 
trary but motivates from an iron calorimeter detector 
planned in India [25| . However, we would also like to 
emphasize that the features of the results and conclu- 
sions obtained below are going to remain the same for 
some other choice as well. Furthermore, in the study 
below, we do not commit to very specific topology of the 
detector and the surroundings. Therefore, one can ex- 
pect some calculable changes in the results for specific 
details of the detector considered. 

It is noteworthy that the underwater or under ice ex- 
periments like BAIKAL [26] and AMANDA are capable 
of, and already measuring the zenith angle and depth 
dependence of the (integrated) muon flux. This is to be 
contrasted with the proposal to employ pair meter tech- 
nique to measure the energy spectrum of muons, with the 
latter giving direct information about energy of individ- 
ual muons. Measurement of muon energy spectrum will 
be of immense use in directly addressing issues related 
to cosmic rays and also high energy neutrino physics. 
We believe that a combination of results from both type 
of methods will play a decisive role in settling various 
issues. 

Muon fluxes in different models of charm pro- 
duction As we discuss in the introduction there exist 
a host of studies on the prompt muon fluxes Since it is 
not possible to provide an exhaustive comparison in this 
letter, we choose to survey some of the cases. 

Following Q the conventional flux (CONV) and the 
prompt flux (TIG) have been parameterized as 



dN 
~dE 



n q e 



E<E a 

1+ AE- E > E a 



1+AE 

n'e-~i 



(4) 



For the conventional muon flux Nq — 0.2, N Q — 0.2, 7 = 
1.74, 7' = 2.1, E a = 5.3 x 10 5 GeV, A = 0.007 while for 
the prompt muon flux N a = 1.4 x 10~ 5 , N' = 4.3 x 10~ 4 , 
7 = 1.77, 7' = 2.01, E a = 9.2 x 10 5 GeV, A = 2.8 x 10~ 8 . 

Gelmini, Gondolo and Varieschi (GGV)[llJ have in- 
cluded NLO corrections for the charm production with 



xg(x) 



, (A varying in the range — 0.5). Their re- 



sults obey th following parameterization for the sea level 
muon fluxes: 



M E ») = A 



E„ 



- [a-\-bx-\-cx -\-dx ) 



(5) 



1 GeV 

where x = Log w (E M / 'GeV) . The parameters are given 
in the Table 1. We choose two representative sets cor- 



responding to A = 0.1 (GGV1) and A = 0.5 (GGV5). 
As our last representative model, we consider flux calcu- 



Flux model 


A 


A xhT 


a 


b 


cxl0~ 2 


dxl0~ ; ' 


GGV1 


0.1 


3.12 


2.70 


-0.095 


1.49 


-0.2148 


GGV5 


0.5 


0.58 


1.84 


0.257 


-4.05 


2.455 



TABLE I: GGV parameters for the prompt muon and an- 
tiunion fluxes. MRSL PDFs are adopted and the scale is set 
to M = 2/t = 2E T = 2(fc|, + ml); m c = 1.25 GeV and k T is 
the transverse momentum of charm quark. 

lation within the saturation model proposed by Golec- 
Biernat and Wuthsoff [U? ]. We follow [12] m parame- 
terizing the charm production within this model. Within 
10% accuracy, we can adopt the following parametriza- 
tion 

x^-Op + air^c+ ..) = AxP{l - x 1 - 2 ) 11 (6) 
ax 

where f3 = 0.05 - ln(£yi0 TeV) and 

0.025 ln(.E/10 TeV) 10 4 < E < 10 8 GeV 



7.6 
7.6 



0.012 ln(£;/10 TeV) 10 8 < E < 10 11 GeV 



(7) 



A = 



140+ (llln(£;/0.1 TeV)) 
4100 + 245 ln(£yi0 5 TeV) 



1.65 



10 4 < E < 10 8 GeV . , 
18 < E < 10 11 GeV W 



10 8 



A is expressed in nano-barns. For this model, we con- 
sider two cases: GBW1 - where the protons are taken 
to be the primary and GBW2 - where we include the 
effect of heavy elements also. On the basis of arguments 
presented earlier, these two cases should have a different 
nature, with the expectation that GBW2 should have a 
decreased number of muons. 

Figure 1 shows the results for muon fluxes in various 
models. In Figure 2, we present the number of muons 
entering a typical detector per solid angle in a span of 
five years. To have a quick comparison, we present the 
same results in Table 2. 
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FIG. 1: E„ 3 x Flux vs. muon enegy in different flux mod- 
els entering the underground detector after passing through 
the rock depth of 3.5 x 10 5 gm/cm 2 . 

Discussion and Conclusions In Table 2 (and 
Figure 2), the number of muons entering the detector 
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FIG. 2: Number of muons entering the 50 kT detector in 5 
years per solid angle vs. for different flux models. Energy 
loss in the surrounding rock of depth 3.5 x 10 5 gm/cm 2 is 
taken into account. 
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TABLE II: Number of muons per solid angle entering the 
detector in 5 years for various energies of the entering muon, 
E„ (in TeV). 

are quoted as number per steradian. However, a bet- 
ter quantity would be the total number of muons. Very 
naively, for the down going muons that we are inter- 
ested in here, this would require multiplying these num- 
bers by a factor 2tt. However, in reality this would per- 
haps be overestimating the number of muons. Therefore, 
a conservative multiplicative factor can be #7r where 
1 < # < 2. From the figures and the table, it is clear 
that direct measurement of muon spectrum using the 
pair meter technique is capable of distinguishing between 
various models of charm production. In particular, as 
expected, making the composition heavier reduces the 
number of prompt muons (the curve GBW2 in Fig 2). 



The conventional contribution is dominant in this re- 
gion. However, due to the fact that the conventional 
flux is much better known and the theoretical uncer- 
tainties arc in better control there, this should be not 
a serious trouble. As far as the prompt muons are con- 
cerned, the conventional one should be viewed as a calcu- 
lable background. Further, from Table 2 we notice that 
Eft ~ 200 TeV seems as a practical cut-off within the 
set up assumed here. However, let us remind ourselves 
that E™ rface ~ 5£ M and E CR ~ 20E™ r f ace and the 
fact that very high energy muons do not lose significant 
amount of energy while traveling downwards. There- 
fore, a muon of given energy entering the detector is 
actually probing the parent cosmic ray of energy almost 
two orders of magnitude higher. We can thus expect 
to probe the hadron dynamics in the region around the 
knee, and also the possible change in composition, very 
easily with this technique. Before concluding we would 
like to mention that there is an indication of strangeness 
enhancement in high energy heavy ion collisions from 
RHIC and CERN-SPS (for a survey of results, see [13). 
Taking clue from this, we expect the contribution from 
the strange hadrons to increase. In particular, for the 
heavy primary composition, the kaon component of the 
conventional flux and prompt component from D s is ex- 
pected to change, though a detailed calculation is needed 
to quantify it. 

In conclusion, we have explored the possibility of us- 
ing pair meter technique to measure the prompt muon 
spectrum. This can be used to gain information about 
the charm production mechanism at these energies and 
therefore yield invaluable information about low-x be- 
haviour of gluon distributions. The range probed here is 
very different from what is expected to be probed at the 
LHC and therefore this will provide an important set of 
data points in the gluon distribution evolution plot. The 
results obtained are quite encouraing and indicate that 
prompt muons (leptons in general) have the capability 
of probing the charm (heavy hadron) production models 
and also the hypothesis of composition change around 
the knee region. 
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